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LANDING-GEAR IMPACT 
By W. Fliigge 


SUMMARY 


This report deals with the impact forces in landing gears. Both 
the landing impact and the taxying impact have been considered, but drag 
forces have been so far excluded. The differential equations are devel- 
oped and their numerical integration is shown, considering the nonlinear 
properties of the oleo shock strut. A way is shown how the dimensions 
of the metering pin may be determined from a given load-time diagram. 

A review of German literature on landing-gear impact is also presented. 


INTRODUCTION 


The objective of this report is to study the impact forces acting 
on the wheels and shock struts of an airplane. For practical reasons 
the investigation has been limited to the vertical forces and does, not 
consider the effect of the drag load which, acts on the wheels during ■■ 
the spin-up time. Within the limits drawn by this restriction, an attempt 
has been made to develop a method for numerical computations which, it 
is hoped, will be useful in practical design work. 

The oleo -pneumatic shock strut which is now in general use and 
which has attained a high degree. of perfection exhibits a rather com- 
plicated relation between the force, the stroke, and the rate of stroke. 
For practical work, it is imperative to express this relation in mathe- 
matical form and to develop a method for the numerical solution of the 
ensuing differential equations.. A detailed discussion of this subject 
will be found in the section "Intentional Nonlinearities." 

Nevertheless it is sometimes useful to consider a highly idealized 
type of landing gear, which has linear differential equations-. Although 
such a model will never correctly reproduce the details of the real 
landing impact, it admits of easy mathematical treatment and permits 
study of questions of a more general character. This has been done in 
the section "Linear Spring -Damper Systems" and the usefulness of the 
results obtained there lies in the fact that they do not depend on the 
more or less incidental details of real landing gears which unavoidably 
enter the computations of the nonlinear theory. 
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The later sections of the report are devoted to refinements of the 
method. Except for the simple mass corrections., they will not be used 
in daily routine work, but they may be of particular importance when 
airplanes of unusual design are built. 

This, work was conducted at Stanford' University under the sponsorship 
and with the financial assistance of the National Advisory Committee for 
Aeronautics. 

While preparing this report the author has received valuable 
information on current American, practice through Mr. J. F. McBrearty, 
Lockheed Aircraft Carp., Mr. K. E. Van Every, Douglas Aircraft Co., Inc., 
and Mr. A. I. Sibila, Chance Vought Aircraft, for which he wishes to 
express his thanks. He also wishes to thank Mr. C. W. Coale for his 
active help throughout the preparation and the writing of this report. 


SYMBOIS 


A n , B n coefficients 

Ao inner cross section of barrel at oil level 

Ap total cross section of piston 

A2 inner cross- section of piston 

A3 area of gap between metering pin and edge of orifice 

a, b distances of landing gears from center of gravity (used 

only in section "The Airplane as a Whole") 

b damping constant for one shock strut (used only in section 

"Linear Spring -Damper Systems") 

F 0 . force in strut when strut is fully expanded and at rest . 

Fp compressive force in shock strut 

F 2 compressive force between wheel and ground (if different 

from Fp) 

F3 force in auxiliary landing gear 

g acceleration due to gravity 
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*2 

m 

m l 

nig 

Po 

Pi 

P 2 

T 

t 

V 

v 

W 

W 1 

x 

X 1 

X 2 
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height of obstacle encountered during taxying 

moment of inertia of airplane with respect to longitudinal 
and lateral. axis, respectively 

radii of gyration 

spring constant for one shock strut 
spring constant for one tire 
mass of airplane 

that part of mass m attributed to one landing gear 
unsprung mass for one landing gear 

pressure in both chambers when strut is fully expanded 
and at rest 

pressure in upper chamber of strut 
pressure in lower chamber of strut 
reference time 
time 

vertical velocity of landing gear when it first touches 
ground 

oil velocity in orifice 

vertical force, other than impact force, acting on 
airplane (weight minus lift) 

that part of force W attributed to one landing gear 

stroke of shock strut 

vertical displacement of mass mp 

vertical displacement of unsprung mass m 2 
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k 

x^ used instead of x-j_ when a second landing gear must he 

considered simultaneously 

z displacement of center of gravity of airplane 

a angular displacement (angle of pitch) of airplane in its 

plane of symmetry 

7 ratio of specific heats 

8 static deflection of mass 

p density of oil 

All displacements are zero when the wheels touch the ground without 
pressure . 


LINEAR SIRING- DAMPER SYSTEMS 
Differential Equations 


Essentially, a landing gear consists of a shock strut and a wheel 
with a tire. The shock strut may he compressed considerably. It opposes 
this deformation with an elastic force increasing with increasing stroke 
and with a damping force which depends on the rate of stroke and which 
dissipates mechanical energy. This shock strut may be represented by a 
spring and a damper arranged in parallel (fig. 1). The tire is for the 
present purposes a simple spring whose deformation is more or less pro- 
portional to the applied force. _ 

Between these two deformable elements there is the mass mg of • 

the wheel, including those parts of the shock strut which participate in 
the motion of the wheel. On top of the whole landing gear there is the 
airplane mass or, more exactly, that portion m-j_ of the airplane mass 

which belongs to the landing gear under consideration. 

When the airplane lands, this system approaches the ground with a 
considerable velocity. As long as the spinning up of the wheels is not 
considered, only the vertical component V of this velocity is of 
interest. The impact begins when the lower end of the landing gear 
touches the ground. This instant is designated t = 0, and the vibra- 
tions of the masses mj_ and mg are studied’ which follow for t > 0 
when the motion of the lower end of the spring-mass system is suddenly 
stopped. 



Linear differential equations are obtained when the shock strut is 
replaced by a simple spring and a viscous damper (dashpot) and when 
the tire is assumed to be a linear spring. There might also be a linear 
damper coupled with the tire but, compared with the shock-strut damping, 
the contribution of the tire to damping is so small that it does not 
seem worth while to include it in the equations. 

The differential equations of the landing-gear impact will now be 
formulated. Figure 2 shows the mechanical system in two positions, one 
for t = 0 and the other for some later time. The displacements of 
the masses m^ and mg, measured from their positions at t = 0, are 
called xj_ and Xg, respectively. Their difference 


x 


= X 1 “ *2 


is the stroke of the shock strut. 

If kg is the spring constant of the tire, then the force trans- 
mitted from the ground to the unsprung mass mg is 


Fg = kgXg (la) 

On the other hand, the force in the shock strut is the sum of an elastic 
force k^x and of a damping force which, in linear theory, must be assumed 

proportional to the velocity x = dx/dt with which the masses m^ and mg 
approach each other: 

F-j_ = k^_x + bx (lb) 

The third force is the load which acts as an external force .on the 

mass m^. It is a part of the weight of the airplane minus a corre- 
sponding part of the wing lift. It will' be shown in the section "The 
Airplane as a Whole" what part of the total weight and lift must be 
attributed to each landing gear. 

The three forces Wj_, F]_, and Fg determine the motion of the 
masses m^ and mg, according to the equations 
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When F-j_ and F 2 are expressed here by x and X 2 according to 
equations (la) and (lb), or better still by x^ and Xg, the differ- 
ential equations of- the linear landing gear are obtained: 


mi’xi + b^Xp - ig) + k-^x^ - x g ) = W X 
nigXg - b^ - Xg) - h- L ^x 1 - x 2 j + = 0 


(2a) 

(2b) 


The problem is of the fourth order and requires four initial conditions 
for t = 0. Before its solution is given, a simplified version will be 
considered which is sufficient in most cases. 


Solution Neglecting Unsprung Mass 


The unsprung mass mg is rather small, usually between 2 and 
5 percent of the mass m^. Under certain conditions it has a very 
definite influence on the force in. the landing gear. But it will be 
seen that it is only of minor importance for the early phase of the 
landing impact, up to and beyond the maximum of the impact force. One 
may therefore begin with a simplified set of .differential equations, 
obtained from equations (2) by dropping the term with m 2 : 


mp'xp + b^X]_ 
-b(xi - x 2 ) 


* 2 ) + k l( x l - x 2) = W 1 

(3a) 

k-^x-^ - x 2 ) + kgXg = 0 

(3b) 


Since x-[_ and x 2 are the highest derivatives occurring in these 

equations, the problem is of the third order, and there must be three 
initial conditions . . 

Two of them follow from the, fact that the displacements xp and x^ 
are counted from the position of the system at t = 0. Therefore 


t = 0: 


xi = 0, 


x 2 = 0 
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The third condition is that the mass has at this time the veloc- 

ity V: 

t = 0: x 2 = V 


It is useful to replace the second condition by an equation for x^_. 
This may easily be done by introducing all three conditions in the dif- 
ferential equations. Equation (3b) yields xg = V, and it is seen here 
that this is not an independent fourth condition, as one might feel 
inclined to think. Equation (3a) yields now: 


t = 0: x’i = W^/m^ 

and this initial relation may be used instead of any one of the other 
three, preferably instead of X 2 = 0. 

One may easily find a particular solution of differential 
equations (3): 


X 1 


= Wl 


k l + *2 
*1*2 





w 


It describes the position in which the system is in equilibrium under 
the load W]_. Besides this, the solution of the homogeneous equations 
is needed. Since all coefficients of the equations are constant, the 
homogeneous solutions are exponential functions of time, say 


X 1 


Ae‘ 


\t 


Xg = Be kt 

When this solution is introduced into equations (3a) and (3b) after 
dropping there the term W-]_, two linear equations are found for A 

and B: 
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. + b\ + k^ - B^bk + k-^ = 


A(b\ + kj - b(' 


"bA, + kn + 


^2) = 


(5) 


Since these equations are homogeneous, they will not have a solution 
A /= 0, B / 0 unless the determinant of the coefficients vanishes, and 
this condition yields the characteristic equation of the problem: 


x 3 + k i + H + h x + 


X c + 


m-, 


m-]_b 


= 0 


( 6 ) 


It is of the third degree. One of its three roots must be real and, 
since all coefficients are positive, this root must necessarily be 
negative, say 


\ - -I3 

The other two “roots may also be real and negative, 

<r 

X ■“ — X-j 


and 

X “ 

or they may be conjugate complex: 


X = -p ± iv 

It may easily be shown that in this case the real part must be negative, 

- 14 . 


For the rest of the formal treatment the cases of real and of „ 
complex roots X must be separated. 
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High, damping, all roots real .- In the case of high damping with 
all roots real the general solution of differential equations (3) is 


k n + 


X 1 = 


^2 —In t -Xot — Aqt 

• A - 1 ■ A 2 e * + A 3 e J 


i i-L-i 

k !^2 


+ A-,e + 


X 2 = ^ + Bpe + B^e + B^e 3 


(7) 


In these formulas only the constants Ap, Ag, and A^ may he chosen 
arbitrarily, while Bp, Bg* and B^ depend on them through equa- 
tions (5), in which in each case the appropriate \ must be inserted". 

When t is set equal to zero and then Xp and its derivatives are 

introduced in the initial conditions, a set of three linear equations is 
obtained for Ap, A £ } and A^. They are: 

A i + A 2 + a 3 = " w l^ + ^ 


IpAp + X-2^-2 ^3 a 3 = 

X 1 2a 1 + ^2 2a 2 + X 3 2a 3 = 


They must be solved numerically, and 
found for any time t. 

The most interesting quantities 
the impact force Fp. For both the 


then, the displacements may be 

are the stroke x = Xp - X 2 and 
B's' are needed in terms of the A's 


Bn = A. 


kp - bl n 
n k-^ + 1<2 - bX n 


n = 1, 2, 3 (8) 


and then 

x . h + Ms + ^ e -V + 

k-L k^ + k2 — "b^.2. + kg — 1 dA, 2 k^ kg "" "^^3 


Y 
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and 


*1 = F 2 


= kgXg 

( ~X.lt -X.pt ~X.pt | 

B 1 e 1 + B 2 e d + B 3 e ^ J 


The coefficients B n need not necessarily be positive. Then there 
may occur a time t > 0, where F 2 = °* "t* 13 - 8 happens, it would 

terminate the domain of validity of the formulas. For greater values 
of t, the force F 2 would not become negative (i.e., tensile), but 

the wheel would leave the ground) the airplane would rebound. 

Because of the force Wi, the airplane would soon return to the 

ground. Meanwhile its horizontal speed or the angle of attack might 
have decreased and hence Wp increased. The vertical velocity at the 
second impact would be, on the other hand, considerably smaller than V. 
The new impact would therefore be less violent, but not necessarily 
uninteresting, because it would find the shock Btrut in a lesB favorable 
condition, with x >0 and, perhaps, close to the possible limit. 

Whether rebounding will occur and how strong the second impact will 
be can be determined only from detailed numerical computations in each 
particular case. But one may say quite generally that the probability 
of a zero of F 2 ,_is greater the more solutions the homogeneous equations 
have - Since each additional mass and each additional spring increases 
the order of the equations, one should avoid mechanical complexity if 
rebounding is undesirable. 


Low damping, one pair of roots complex .- The complex exponentials 
which appear in the case of low damping with one pair of roots complex 
may be expressed in real form by exponential and trigonometric functions: 


= W, 


ki + 


kg 


kikg 


+ e 


-M-t 


{ A 1 


cos vt + A 0 sin 


vt) 




* W 1 4 + 


cos Vt + B 


2 


_ain vt) 


+ B^e 


-X.pt 


( 9 ) 
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The relations between the A's and the B's are here more involved 
because equation (8) can be applied only before the trigonometric 
functions are introduced. When this is done the following relations 
are found: 


Bp = (l + a i) A l + Pi^ 
B£ = - Pi A i + (l + a 2_)A2 
B^ = ^1 ^ ■^‘3 


with 


a ! = m l 


P 2 ^ - bp) - v 2 ^ + bp) 
(k]_ - bp) 2 + b 2 !' 2 


Pi =m 1 


bv(p 2 + v 2 ) - 2kjPV 
^k-j_ - bp) 2 + b 2 v 2 


m-, 1 


Og = 


13 


k-^ — bl 3 


The boundary conditions will now yield the following set of three 
equations for A^ Ag, and A^: 


A-i + Aq = -W if— + — ) 

1 3 H*i *e 


pAj_ - VAg + X^ a 3 = “V 


(p 2 - V 2 )a! - &1VA 2 + X 3 2 A 3 = ^ 
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When these equations have been solved numerically, B]_, B 2 .> 

and B3 may ,be found from the preceding formulas. Then 


X = X^ - Xg 


% _ A 3 m l x 3 2 e “ x 3 t 


(“tAl + ^r^) 6 ^ 


COS Vt + 


( P 1 A 1 “ “ l ^) 6 tit 


sin Vt 


( 10 ) 


F, = Fr 


— h^x^ 



+ 




cos vt + B 2 


sin vt 


l 


In this case rebounding is rather probable because of the trigonometric 
terms, and it will be inevitable when = 0. 


Solution Wot Neglecting Unsprung Mass 

The results just described may be considered representative for the 
landing impact if it can be shown that they are not seriously affected 
by the neglected mass m2 of the wheel. This- side of the problem will 
now be investigated. 

Instead, of equations (3) set (2) must be used which still contains 
the term with mg. Four initial conditions are required: 


t = 0: x-j_ = 0, x 2 = 0, Xj_ = V, x 2 = V 
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which now are all independent of each other. When they are introduced 
into differential equations ( 2 ), there are obtained 


t = 0: Xi = Xg = 0 


If equation (2a) is differentiated once more and again everything known 
is introduced, another dependent condition is found: 


bW, 

t = 0: % = - — I 


For actual use, choose from all these conditions the set of four which 
refer to x-j_ only: 


W, 


t = 0: x x = 0, x ± = V, x x = — , 


bW, 

_ 1 

12 

m l 


( 11 ) 


The general solution of equations (2) is 


k x + k2 -p t N 

1 = W 1 — k ' kg" ' + ® (Ai cos Vl t + Ag sin v-|_tj + 


-M-ot . V 

e ^ ( a 3 C0S v 2 ^ + a 4 sin v 2 ' t ) 


(12a) 


1 “M-nt, N 

Xg = Wf + e 1 ^ cos V x t + Bg sin V^tJ + 


cos V gt + sin Vgt) 


(12b) 
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Here p and v are obtained from the real arid imaginary parts of the 
four solutions 


X, = -4 2 ± iv-j_ 
X = -P 2 ± iv 2 


of the frequency equation 



f 1 

+ 1\ 

13 / 

+ ( 


+ m 2) 

jbX + I 


m n 


mg 


U2 + l^. Xt l*,o (13) 
/ num, 


■ 2^2 m l m 2 


and the B 1 s depend on the A's by the relations 


with 


Bn = 


= 


B„ = 


3 

B4 


^1 + °-]_)aj + PpA 2 
-Pi A i + (l + a 1 )A 2 

^1 ■+■ ^2) ^ 2^4 

* P 2 A 3 + i 1 + a 2) A 4 




a n = m l 


L n{\ - bp n ) - Vn 2 ^ +_bp n ) 
( k l “ ^n ) 2 + ^n'~ 


Pn = “l 


bv n(' I n g + v n g ) ~ gk llVn 

(*1 - b ^n) 


2 - b2v n 2 


n = 1,2 


(14) 


( 15 ) 



All these formulas are established In the same way as are the corresponding formulas for mg = 0. 

It may happen that one pair of X's is real (or even both pairs). The formulas will then 
undergo similar changes, like those explained for the third-order problem. 


When solution (12) 1 b introduced into the Initial conditions (equations (ll)) the following 
equations for to will result: 



~^1 A 1 + v 1 a 2 “ |I 2 A 3 + v 2 a 4 = v 

(^l 2 - Vi 2 )a! - S^Ag + (n 2 2 - V 2 2 )a 3 - 2^ 2 V 2 A4 = ~ 

hi^i 2 - 3 V! 2 ) A! - - V! 2 )a 2 + h 2 (h 2 2 - 3V 2 2 )a 3 - Vg^ 2 2 



(16) 


This set must be solved numerically and then one may find B-^ to from equations (l^) 

and ( 15) • The stroke x = xp - x 2 of the shock strut Is then: 


Wn 




|( a i A i + ^i^) COB v i^ ” (^i A i ~ BiEL v i t 




(agA 3 + PgA^ c o 0 V 2 t - ^PgA 3 - a^A^ sin Vgt 


(17a) 


H 

vjl 
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The force acting between the tire and the ground is again 


f 2 = ^ 2 X 2 

and may be computed from equation (12b). It is responsible for the 
stresses in the tire, and when it becomes zero the airplane will rebound. 
But F 2 is not equal to the shock-strut force F]_ which in turn is 
responsible for the action of this part of the landing gear and for the 
dynamic load on the airplane structure. This force must be found from 
equation (lb): 


(bUl - kp) ( c^A-l + 3 x A 2 ) + bV 1 (p i A 1 - c^AgM cos V-jt + 

|^ v l( a l A l + Pl A 2) " (^1 “ k l)(Pl A l - a l A 2)] sln v lt j- + 

e - k l) ( a '2^3 + $2^k) + bv 2 (^^3 “ a 2 A 4)J cos v 2 t + 

jbv 2 (a 2 A 3 + 0^) - (bp 2 - k x j (0^3 - a-gA^Jj sin v g t i (17b) 


F x = W x + e 


-p x t J 


From the formulas shown here it is clear that the final result is 
connected with the data of the problem through" an algebraic equation of 
the fourth degree and through a set of four linear equations. Solutions 
must be obtained numerically for a given set of data, so it is not pos- 
sible to discuss the features of the solution in general terms. To find 
out how they look, a series of systematically chosen examples has been 
computed which will be discussed in the section "Discussion of Numerical 
Results." 


Undamped System 

As a basis for this discussion, it is useful to consider the case 
when the damping b is zero and, additionally, mj_ » m 2 . The second 

assumption is certainly good but, if used alone, it would not give any 
substantial mathematical relief. The first assumption is, of course, 
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not very realistic in a landing gear whose essential purpose iB damping, 
but the conclusions derived will help in understanding the more realistic 
cases. 

When b = 0, the frequency equation (equation (13)) loses the 
terms with X^ and \ and, when it is assumed that m^ » mg, the 
first part of the coefficient of X^ may be neglected. 

The equation reads then 


+ fol ia ^ + fofo = 0 


Ulp 




and its solutions are purely imaginary, say 


and 


X = tiVn 


X = ±i Vr 


When the equation for X 2 is solved and everywhere l/mp is neg- 
lected against l/m2, it follows that 


X 


2 


k l + *2 + 


2mg 


+ kg 
2m2 


The upper sign yields X 2 = 0 


and hence v-j_ = 0> the lower sign yields 


v 2 = 


l k l + ^2 

nig 


This indicates that « Vg, but evidently too many small terms were 
neglected to find a reasonable value -for V-i. It may be obtained from 

p p 

the fact that the third term in the equation must equal v ]_ ' 



k l*2 


m^mgVg 


2 
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V -i = 


k l*2 


1 y m l( k l + 


This is indeed much smaller than Vg. 

Since = M-g = Oj solution (12) loses the damping factors e 
and becomes: 


k l + ^2 

X 1 = W 1 ""k kg + A 1 cos v l^ + ^2 sin v l^ + a 3 cos v g' fc + sln v gt 


x 2 .= ^i + Bj_ cos v^t + Bg sin v-j_t + B^ cos Vgt + B^ sin Vgt 


When this expression for x-]_ is introduced in initial conditions (ll), 
a very simple set of equations is obtained: 


Ai + A.3 = 




V-jAg + VgAj^ = V 


5 P «1 

v 1 2a 1 + v 2 2a 3 = - 51 

v 1 3 A 2 + V^A^ = 0 


> 


(18) 


It consists of tvo independent pairs which will be solved and discussed 
separately . 
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From the first and third equations: 


A 


1 



*l + * 2 \ 
* 1*2 ) 


m 2 

n m l (*1 + *2) 2 

k 1 + 3^2 m]_ 

M 1 

IS * 

1 

H 1 


« -W-l 


k-L + 
* 1 * 2 _ 


A 


3 




*1 + * 2 \ 

* 1*2 ) 


W, m 0 

* " E ^ 1 " 1} 

= 0 

Equations (15) yield P-j. = P 2 = 0 an( l 


and equations (lk), 





= -Wi/k2 





= 0 


n = 1, 2 
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The last equation is not very convincing as it stands, since it gives 
as a product of Ay which is almost zero, and of the large 

factor mi/nig. But one may check the result by first finding an exact 
expression for B 3 and then neglecting .against m^. 

From the results it is seen that only the low-frequency motion 
(Ai, B^ is of importance and that the ratio x^: x^ is at all times 

the same as for the static deflections, uninfluenced by the presence 
of the unsprung mass. 

The second and fourth of equations (l 8 ) may be handled in the same 
way. The result is this: 



B 4 " 



= V 


k^ + kg 


Here again it is seen that for x-^ the low- frequency motion 
is by far preponderant, A^ being smaller by a factor 
but for the displacement of the wheel ^Bg> the factor is only 
(mg/mi) 1 / 2 . In the low-frequency motion the ratio xg: Xp is the same 
as for the static deflection, but in the high-frequency mption x g is 
much larger than Xp the wheel moving up and down between the ground 
and the almost unmovable airplane mass. 


On the whole this analysis shows that through the presence of an 
"unsprung mass" a high-frequency motion is added to the low-frequency 
motion of the airplane. This high-frequency motion does not affect x 


1 ' 
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"but it makes a certain contribution to the forces in the shock strut and 
in the tire. Since ^mg/m}) 1 / 2 is still as small as 0.15* one may- 
neglect the presence of the mass mg if the accuracy- requirements are 
not too high. In view of the many arbitrary assumptions which enter the 
analysis (e.g., the value of V) , one may think of neglecting m 2 for 

design purposes, but one should keep an eye on it when evaluating tests. 
This rule, of course, is derived from the behavior of the undamped shock 
strut. How far it is modified by the damping can be seen only from the 
systematic numerical work which will now be discussed. 


Discussion of Numerical Results 

Dimensionless parameters . - The formulas developed in the preceding 
sections have been used to compute some typical examples. In order to 
draw maximum information from this work, it has been done in dimension- 
less form, and therefore the choice made for the dimensionless quantities 
must be discussed before the results may be discussed. 

For the displacements x-^ and xg and the stroke x a reference 

length is needed, and when they are plotted against time a reference 
time is needed. Since the deflections start from zero at t = 0 and 
approach asymptotically definite values, the static deflections, it 
seems reasonable to adopt the static deflection of the mass m-^ aB a 
standard of length: 


8 ' mi % + k) 

A simple time standard may be found in the period of the vibrations which 
the mass m^ can make on the springs k-j_ and kg in the absence of 

damping. This period is 2jt/s/g} drop the factor 2n and choose as 
time standard 



' m l(^i + kg) 

*1*2 


For the forces F^ and Fg the load W-j_ might be used as standard} 
but, since W-]_ depends on the horizontal speed of the airplane and on 

the angle of attack of its wings, it may have rather different values for 
different landing cases of the same airplane. It is therefore better to 
use m-^g as a reference value. 
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Besides these variables there is a set of constants which influence 
the impact. They are the masses m^ and m g, the spring constants k^ 

and kg, the damping b of the shock strut, "the load Wj_ (weight minus 

lift), and the vertical velocity V of the airplane. These constants 
may be combined to the following dimensionless parameters: 

£ — bT/m^, W^/m^g, and VT/S. These parameters must be chosen 

for each example . 


Influence of unsprung mass .- Of most interest in a study of the 
linear spring-damper system is the influence b)f the unsprung mass on 
the impact force. Since this influence is small in the undamped system, 
one may hope to find the same result when damping is present. To check 
whether this is true, a landing gear has been investigated analytically 
for two extreme values of- the mass parameter, nig/mp = 0 and 
1112 /m 1 = 0.050. Values common in current practice lie approximately 
halfway between, and the choice has been made in order to make the effects 
more clearly visible. For the other parameters the following values were 
chosen: 


k l 

k-|_ + kg 


0.25 


bT/m-j_ = 0.5 


Wi/mig =0.2 


VT/S =2.0 


The weight parameter lies halfway between a fully buffered landing with 


= 0 and the usual assumption of 


Wn = — 


m^g. The other three figures 


are so chosen that they correspond to an actual airplane, at least so 
far as a correspondence between a linear and a real shock strut is 
possible. 


The result of the computations is seen in figure 3(a) which shows the 
shock-strut force_ Fp against time t. At the start there is a definite 

difference: The unsprung mass absorbs the first impact. and the shock- 

strut force develops slowly; a slight overshoot follows; and then the 
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curves are practically parallel, the distance between them corresponding 
to the addition of the wheel mass mg to the total mass m^ + mg which 

must be decelerated. 

This result confirms the view that the maximum of the shock- strut 
force is not much influenced by the unsprung mass and that one may 
safely assume mg - 0 when this leads to a simplification of the 
theoretical or numerical work. However, this simplification may not be 
admissible when the wheel spin-up must be considered. Figure 3 ("b) shows 
the force Fg between the wheel and the ground, and here it appears 

that in the early stage there is a considerable difference between the 
two cases such that Fg and hence the drag load will increase with the 

unsprung mass. 

Influence of damping .- High damping in the shock strut is desirable 
since it dissipates the kinetic energy of the airplane and thus prevents 
repeated rebounding. The influence of damping on the landing impact may 
be seen in figure 4. Here the force Fg on the wheel and the stroke x 
of the shock strut are plotted against time for the following set of 
parameters: 


mg/mi = 0.025 


kf + kg 


0.25 


w l/ m lS = 0.2 


VT/5 = 2.0 


and bT/m-L = 0*5 and 1.0. In the initial stage there is not much dif- 
ference between the two curves (fig. 4(a)), since the impact is caught 
by the tire, but then the rise of Fg is much faster for the case of 

higher damping; the maximum is reached more quickly but it is only 
7 percent higher than that for the case with half as much damping. The 
development of the stroke (fig. 4(b)) is on the whole similar in both 
cases, but the maximum is lower for high damping. The figures show 
that, apart from its influence on rebounding (which occurs much later), 
high damping has its pros and cons and that they must be balanced in 
each design. 
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Influence of spring constants . - In another example the influence of 
the spring' constants k^ and kg is compared. The forces F^ and Fg 

have been computed - for 

m 2/ m l = 0-025 
bT/mp = 0.5 ~ ' 



VT/6 =2.0 


and for — — = = 0.15 and 0.25. ’ 

*1 + k 2 . : 

The standards of length and time, 5 and T, will not be changed 
if the sum of the reciprocals 1/k-, + 1/kg is kept constant, and under 

ki + kg) simply means 

that the shock strut is made stiff er and the tire softer. Consequently, 
a greater part of the total deformation will take place in the tire, and 
since there is no damping the impact force will build up more slowly. 

This is clearly seen in figures 5(a) and . 5("b) , which show that the 
influence is the strongest on Fg. Of course, the smaller the impact 

force, the less the airplane will be decelerated and the higher the force 
must-rise at a later, stage to bring the vertical motion to a stop. This 
is also seen in the diagrams, and in figure 5(c) one sees the consequences 
for the displacement xp of the airplane and the deformation x>> of 
the tire. _ 

Influence of weight and lift . - Remember jthat the notation Wp 

represents the resultant static load on the landing gear, essentially 
the difference of the weight of the airplane and the wing lift. In the 
landing-gear literature one finds discussions of the whole gamut of - 
possibilities from the buffered landing with W-j_ = 0 to the pancake 

landing Wp = mpg. The present American regulations consider 
W/mpg = 0.333 as a standard- assumption. The, cases where Wp/mpg = 0 
and 0.2 have' been computed, assuming :.r 


this condition an increase of the parameter kp/ ( 
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m2/mi 


= 0 


k-^ + kg 


= 0.25 


bT/mp =0.5 


VT/8 =2.0 


The corresponding forces Fg = Fp are shown in figure 6, and one 
recognizes that the difference is not very pronounced. This is easily 
understood when one considers the maximum displacement Xp of the air- 
plane mass. It is 1.655 in the first case and 1.826 in the second, 
while the final static deflections will he 0 and 0.26, respectively. 

It is therefore essentially the kinetic energy of the airplane mass and 
not the weight that is responsible for the impact. 

All the examples given here show the general trend of changes which 
a change of one of the parameters will induce. In the details much will 
he different when real shock struts with their essential nonlinearity 
are considered. 


Taxying 

When an airplane taxies on the ground, Its tires and the shock 
struts have the same functions as the tires and springs of an automobile. 
Whether they will he subject to serious dynamic forces depends on the 
smoothness of the ground. In the investigation of the taxying impact 
it has become customary to assume that the airplane rolls at moderate 
speed over a bump shaped after a sine curve. On a turf -covered airfield 
such a bump may represent a frozen molehill or a similar obstacle, but 
on a well-kept concrete runway it Is difficult to discover an obstacle 
of this kind from which the length and height of the bump might be 
taken, and the same is time for the deck of an aircraft carrier. 

It is preferred therefore to assume as a standard obstacle a step 
in the ground, as it is encountered In the joints between the runway 
slabs or if a wheel should get over the edge of the pavement (fig. 7) • 
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For the mathematical formulation of the. problem it is assumed that 
the airplane is taxying on the left part of the pavement and that the 
landing gear has the static deflection: 


X 1 = W 1 


kj_ + k£> 
k l k 2 


x 2 



Since taxying is done at low speed, in these formulas comes close 

to the total weight m-^g to be carried by one shock strut. 

At t = 0 the wheel hits the step in the pavement, and for t > 0 
the term kpx 2 in equation (2b), which represents the force in the tire, 
must be replaced by kg(x 2 + h) . The equilibrium is then disturbed, and it 
is desired to know the resulting vibration. It will be found by solving 
the differential equations - 


m l x l + b ( x l " * 2 ) + k l( x l ~ * 2 ) = 


Wn 


(19a) 


0 ^X 2 - b - x^ - k-j^Xj - x 2 ) + kr^ = -kgh 


(19b) 


for the following set of initial conditions: 


k, + k 

* - 0: x i = w i 


x 2 " "iX? 


X 1 = *2 = 0 


It will again be useful to write all initial Icondit ions in terms of one 
variable. Since it will be seen that In the present ..case x 2 is more 
important than x^, x 2 is choBen to formulate the conditions. The 

procedure Is almost the same as that described in the section "Solution 
Not Neglecting Unsprung Mass," and the result is thiq: 


k2h 
--mg } 



( 20 ) 


t = 0: 


x 2 - W 1 


x 2 = 0, 


x 2 = 


x 2 = 



4 - 


The solution of equations (19) has, of course, almost the form of equations (12). Because 
of the new term on the right-hand side of equation (19b) it is: 


X 1 


W-, 


4* k 2 


*1*2 


h + e 


-^t 


cos V x t + Ag sin VjtJ + e" M '2 t ^ cos Vgt + sin V 2 tj (21a) 


*2 


= W, 


*2 


h + e 


Bj cos Vit + B 2 sin V 


x t ^ + e ^/ 03 cob v ^t + B]^ sin Vgtj 


(21b) 


When this is Introduced into the initial conditions, equations (20), the following set Of 
equations for Bi to B^ will result: 


= B 


~ V^Bg ^ — ® 


(Hl2 _ ^2j Bi _ - v 2 2 )b 3 - SUgVgB^ = - — 

^l(^l 2 “ 3y! 2 )Bi - V^ 2 - Vi 2 jB 2 + H 2 (^2 2 - 3V 2 2 )b 3 - V,^ 2 - V^Bj, = 

When these equations have been solved, the constants A must be computed, 
inverting equations (lk): 


r (22) 


bk^ 

__ 

m2 

J 

This is done by 


(l + ai^Bi - P x B 2 

Al - (l + a x f + 


ro 

-4 
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(1 + + Pi 2 


A 


3 


(1 + a g )B 3 - 3^4 

(l + c ^) 2 + p 2 2 


PgBg + H- QOgjB^ 
(1 + og) 2 + p 2 2 


Now and X2 and their derivatives may be calculated as functions 
of time and from them,, the stroke x = x^ - x 2 , the force on the wheel 
F 2 = kg^Xg + h), and the force in the shock strut 

Fj = - m^x*2_ 



An example may illustrate the mechanical content of these formulas. The 
following set of dimensionless data is chosen: 
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“2/®! = O' 025 


± = 0.15 

+ k2 


bT/m-]_ = 0.5 


w l/ m lS = 1 
h/5 = 0.1 


After solving frequency equation (13) one may find B 1 to B^ from 
set (22): 

B ± = 0.01378 
B 2 = 0.0115 

b 3 = 0.08636 

B^ = 0.06176 

The precision of B 2 is rather poor, hut it is not possible to obtain 
a more accurate value unless the data of the problem are given with such 
accuracy that a computation with more than slide-rule accuracy would be 
justified. However, B^ is multiplied in equation (21b) with a factor 

which increases rather slowly with t, and the term does not reach an 
important magnitude before the essential phase of the impact is passed. 

The displacement X2 of the wheel is shown by the solid line in figure 8. 
The displacement x-^ of the airplane has not been plotted because it is 

almost constant, and only after a considerably longer lapse of time does 
the airplane climb Blowly to the new level given by the step in the 
runway. 
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The forces Fp in the shock strut and Fg in the tire are 
represented Toy the solid lines in figure 9- The force Fg jumps _ 

instantaneously from its static value Fg = Wp to 1.6jV!^, the rise 
being determined by‘"the height h of the step and by the stiffness kg 
of the tire. Actually this sudden increase of F2 is smoothened by 
local deformation of the tire. The force Fp which through the shock 

strut acts- on the airplane structure rises smoothly to a maximum and 
then returns in damped oscillations to its static value Wp. In the 
present example 'the essential part of the impact is passed at t = 0.15T. 

The fact that the mass mp hardly moves, within that time which is 
of interest suggests simplifying the computations by putting mp = °°. 

When this 1 b done, equation (19a). must be dropped entirely (it simply 
yields xp = 0) and in equation (19b) there must be put 


= 5 = W n 


k l + 


2 - 




« 

X 


1 


= 0 


The problem is then reduced from the fourth to the second order and its 
solution is 


W 1 

X 2 " k^ 


kgh 


kp + kg 


+ e 


-Pt 


^Bp cos Vt + B 2 sin vt^ 


with 


P 


b 

2mg 


v = aj fFi - - ” 2 

Of the initial conditions (equations (20)) only the first two remain 
valid, and from them 
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1 Cgh. 

k-j_ + kg 



With these formulas it is rather easy to compute the deflection x 2 
and the forces 


F 1 = k l( x l ~ x 2) + b (*l “ * 2 ) 

kp + ko 

= W 1 — * 1 * 2 . ~ bx 2 



Fg = kgXg 

However, there is still a difficulty in the dimensionless representation 
of the results. When m 1 is put equal to », the quantity which was 

used as a reference basis for mg* bT, and seems to be lost. But 

this difficulty is only apparent. When m^ is set equal to °° in an 
equation, this does not mean that the mass really is infinite but only 
that the inertia is intentionally overrated. Nevertheless, there is a 
certain weight m-j_g to which the load (weight minus lift) W^_ and the 

weight mgg of the wheel may be referred and which will produce a 
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certain static deflection 8. To this deflection the height h of- the 
step is referred,, and the damping h is handled by writing 


bT _ _b_ 

m-j_ m-|_ 



b 



With this interpretation of the dimensionless quantities, the 
results of the simplified theory may be plotted in the same diagrams 
which were used before. They are represented by the broken lines in 
figures 8 and 9« ..... 

These curves show the following features r~ (l) In the domain of 
interest they are so close to the exact curves that they can hardly be 
distinguished. (2) For large values of t they have different 
asymptotes. This is easily explained. When ;x^ = Constant, the springs 

will find themselves at last more compressed than they were before the 
impact. The wheel can therefore not rise by the full height h of the 
step, and the forces F^ = F 2 will be higher than W-j_. However, this 

deviation between the two solutions is of.no practical importance, not 
so much because of its small magnitude but because of~its late occurrence. 

It is of some Interest to study the extreme case that m-^ = 00 
and m 2 = 0. Since the mass m 2 is responsible for the difference 
between the forces F-j_ and r 2; it is seen from a glance at figure 9 

that this simplification of the problem goes too far to yield results 
of immediate practical value. The formulas therefore will not be 
reproduced, but some points computed from them_have been entered in the 
force diagram (fig. 9)* This line of - dots which represent both Fp 

and F 2 shows approximately how the solution will be changed if the 
unsprung mass m 2 is substantially decreased: The sudden rise of F 2 

is the same, but the following decrease is faster. The force F^ rises 
more rapidly ( in the limiting case has the same discontinuous increase 
as Fg).> and its maximum will be higher the more is decreased. This 

shows that in taxying it is not advantageous to have the unsprung mass 
too small. 



NACA TN 274-3 


33 


INTENTIONAL NONLINEARITIES 


The kind of shock strut considered in the section "Linear Spring- 
Damper Systems" is the only one which leads to linear differential 
equations. The spring terms in these equations correspond to the action 
of helical or other steel springs, and such springs were used to some 
extent in early shock struts . Modern shock struts use air as an e la stic 
material^ and air does not show linear elasticity unless there iB time 
enough to dissipate the heat generated by compression. However, the 
nonlinearity introduced by a pneumatic spring is not severe, even in 
the extreme case of adiabatic compression. 

Quite different is the situation with the damping term in equa- 
tion (lb). Viscous damping is never realized in shock struts, their 
damping being produced by the acceleration of oil squeezed through 
narrow orifices or slots. If the cross section of the orifice does not 
vary, one has a velocity-square damping, and this already presents an 
essential nonlinearity. But more than this, the necessity of making 
the best use of the structural weight of the landing gear has led to 
the introduction of a metering pin which changes the width of the 
orifice in such a way as to make the impact force increase quickly to 
its peak value and then stay at this value for a considerable time. 

The nonlinearity which the metering pin introduces into the differential 
equations is intentional and essential, and one has to study the equa- 
tions of motion with the corresponding damping term. This will be done 
in this part of the report and, since there is no additional difficulty 
connected with it, the nonlinear elasticity of the air spring will also 
be included. 


Differential Equation of Oleo Strut 

Oleo struts are built in different forms (fig. 10) . They all have 
this in common: A piston moves in a cylinder, and there are two chambers, 

separated by a diaphragm and connected by an orifice. The lower chamber 
is filled with oilj the upper one, partly with oil and partly with air. 
When the strut is compressed, oil must flow from the lower to the upper 
chamber, and there may or may not be a meter in g pin which fills part of 
the orifice and makes the remaining gap depend on the position of the 
piston. 

The pressure p^ in the upper chamber depends only on the air 
volume and hence on the position of the piston. The pressure p£ in 
the lower chamber is greater by the pressure which is needed to squeeze 
oil through the orifice. The difference is proportional to the square 
of the piston velocity and thus produces a damping of the shock-strut 
motion. 
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The relation between the force transmitted through the shock 

strut, the stroke x, and the rate of stroke x will now be established. 
In addition to the notations explained in figures 10(a) and 10(b), the 
following symbols : for the different cross-sectional areas are used: 


A q -inner cross section of barrel at oil level 

Aj total cross section of piston 

Ag inner cross section of piston in figure 10(b) 

When the strut is fully expanded and at_rest, there will be a 
certain pressure p Q in both chambers and a force up to the limit 


F q = V Q A ± 

may be applied without displacing the piston. 

When the piston is displaced, the content of the chambers is 
decreased by xA^ and, since the oil is incompressible, the air volume 
must decrease by this amount: 


xA x = (z 0 - z)Aq 

The collapse of a shock strut under the landing impact takes less 
than 1 second, and one might think that this time would be too short to 
allow for much heat transfer. Then the compression of the air would 
follow the adiabatic law 


7 7 

Pl z = Po z o 


with 7 = l.k. There is, however, a very efficient cooling of the air 
through the jet of cool oil which is shot vigorously through the orifice 
and scattered on the cylinder walls. It may therefore be justified to 
assume a much lower value for the exponent, say 7 = 1.1 or even 
isothermal compression with 7=1. To decide this point, temperature 
measurements in the air chamber would be needed. 
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Elimination of z from the last two equations yields the following 
expression for the pressure in the upper chamber: 


p l 


= P, 


' z 0 v A i y 

ZqAo/Ai - x) 


( 23 ) 


When the piston is displaced, the volume of the lower chamber is 
decreased and oil must flow through the orifice into the upper chamber. 
In the case of figure 10(a) the rate of the oil flow is Apxj in the 
case of figure 10(b) it is AgX, and the oil velocity in the orifice is 


v 


A^x 

A 3 ~ 


or 


AgX 



Here A 3 represents essentially the area of the gap between the 

metering pin and the edge of the orifice, inclusive of an orifice 
coefficient, if necessary. This gap area depends in a known way on 
the stroke x. But A 3 includes also any other leakage between the 

two chambers, and such additional gaps may depend on elastic deforma- 
tions and hence on the pressure pg. Complications are avoided by 

disregarding this fact and assuming that A 3 Is known or sufficiently 
estimated as a function of x alone. 

The oil is accelerated to the velocity v by the difference 
between the pressures Pg and p^_ in the two chambers according to 

Bernoulli's equation 


P 2 - Pi = fpv 2 


where p is the density (mass per unit volume) of the oil. Taking the 
last two equations together, a relation Is obtained between the pres- 
sures and x: 
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( 24 ) 


respectively, for figures 10(a) and 10(b) . 

In this form, the relation holds only for the upward stroke. 

During the recoil motion of the piston the oil moves in the opposite 
direction through the orifice and the pressure in the upper chamber is 
the higher one. To cover this motion, one must at least write p^ - Pg 

instead of - Pp in equations ( 24 ) . However, even this will not 

really describe the recoil motion, for the following reason: The oil 

jet which is shot in the upper chamber during the upstroke is so 
vigorous that air and oil get thoroughly mixed, and this foam is squeezed 
back during the downstroke. It is therefore scarcely possible to calcu- 
late the details of the downstroke until experimental information becomes 
available concerning the degree of mixing and the density p which should 
be used in equations ( 24 ) for this phase of the motion. 

In figure 10(a) the force acting downward on the piston is simply 
f 1 = a 1 p 2 . 

= A lPl + A x (p 2 - P p) 

In figure 10(b) part of the piston protrudes into the upper chamber 
and is there exposed to the pressure pp. The total force is therefore 

in this case 

F 1 ' ( A 1 • A s) p l + A 2 P £ 



' A l p l + A 2 ( p 2 ' p l) 
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With, equations (23) and (2^) this yields for figure 10(b): 


F 1 = A lP 0 


z o A -o/-^l 
, z o A o/ A l - 




(25) 


and for figure 10(a) the relation is the same except that A 2 must be 
put equal to A]_. Equation (25) is the equation of the oleo -pneumatic 
shock strut. 

For obvious reasons real shock struts differ from the idealized 
figures 10(a) and 10(b) in that the upper end of the piston is so 
shaped that it touches the wall of the barrel. In this way a separate 
annular space is created between the plunger piston and the wall of 
the barrel which is usually connected by good-sized holes with the 
upper chamber. The oil flow through these holes may add some damping. 
It is easily possible to take care of this effect by a correction of 
the factor of the second term in equation (25). 


Dynamic Equations for Landing Impact 

For this study it will be assumed that the force in the tire 
follows a linear law: 


F 2 = kgXg 

but that the force F-j_ in the shock strut depends nonlinearly on the 
stroke x and on the rate of stroke x: 


F 1 = F l^ x ) 


the function F^ being given by equation (25). 

The equations of motion are essentially the same as equations (2) 
except that the terms 


b ^x^ - x 2 ) + k^ (x^ - Xr,) = bx + kpx 
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must be replaced by F-^x, x) . The equations are therefore these: 

m l* x "l + ^]_(x, x ) = ^1 (26a) 

m 2 x 2 _ F i( x .> x ) + ^2 X 2 = 0 (26b) 

Since these equations contain three unknown variables, a third equation 
is needed, the relation 


X = X^ - Xg (26c) 

The initial conditions for these equations will be different from 
those used with equations (2). Because of the prestressing of the oleo 
strut there will be a short but finite time at the beginning of the 
impact when the piston does not move and the total of the deflection 
comes from the tire. During this interval the motion is governed by 
the differential equation 



+ m, 


3) x 2 + kg> 


=-w n 


and to it the initial conditions 


t = 0: 




= V 


must be applied. The solution is 


W 1 V 

x-i = x n = : — (l - cos art) + — sin art 

1 2 ko ao 


kg 


2 


or 



+ mg 
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It is valid up to the time t = t 0 when for the first time 

k 2*2 = F o 


the prestress force of the shock strut. From 



\ V 

- cos ot 0 ) + — sin 0Jt o = F 0 


one finds t and then 


V 

Xq = ^1 - cos 0Jt o ) + - sin 03t c 


Oi 


V 


o 


W^CD 

= -r — sin cnt 0 + V cos cot 
k 2 


o 


the displacement and velocity at the end of the initial interval. With 
these quantities the initial conditions for equations (26) may now be 
written. They are 


x x = x 2 = Xq 

► 


(27) 



J 


These conditions ought to be imposed at t = t 0 . For the numerical 
solution It is more convenient to start a new time scale in which 
equations (27) are to be satisfied at t = 0. When plotting the results, 
one should of course convert the time so that the zero is at the moment 
of first contact. 


Equations (26a) to (26c) with boundary conditions (27) must be 
solved by a step-by-step integration. As is well-known, this is done 
in the following way: At a certain time t, the differential equations 

are used to compute the numerical value of the highest derivative of 
each unknown; then one of various integration methods is used to find 
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values of the unknowns themselves for the time t + At. The choice of 
the integration method determines the exactness of the result and the 
amount of computation work needed. This point will be discussed in 
the next section. 

The other part of the process, the computation of and Xg; 

is done in table 1. 



In the first line of this table, columns ( 2 ) to ( 5 ) are filled 
from the initial conditions. Then columns (6) and ( 7 ) are filled with 
the help of equation (26c) and columns (8) to (10) and (ll) and (12), 
with the help of equations (26b) and (26a). The integration process 
will then yield values for the second line of columns (2) to (5), and 
then the whole procedure may be repeated. ' J. 


Simplified Equations, Neglecting Unsprung Mass 

As has been seen before, the unsprung mass does not essentially 
influence the load-stroke curve. It is therefore of interest to 


reconsider equations (26a) 

to (26c) after dropping the term with 

m 2 : 


m l^l = Wj_ - F|_ 

(28a) 


F X (X, X) = kgXg 

(28b) 


X = Xj_ - X£ 

(28c) 

As seen already in the linear case, neglecting the unsprung mass 
the order of the problem by one and the initial condition for Xg 
be dropped. 

decreases 

must 
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Equations (28) may be handled numerically with the help of table 2, 
combined with a table for the numerical integrations. 


TABLE 2 


(1) 

(2) 

(3) 

w 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

t 

X 1 

x l 

*2 

X 

F 1 

X 

*2 

Wi - F X 

■ • 

X 1 

0 











Here in the first line columns (2) to (4) are filled from the initial 
conditions and the other columns, with the help of equations ( 28 c ) , 
(28b), and (28a). The first step of the numerical integration serves 
to get the second line started, and so forth. 


Methods of Numerical Integration 

There exists much literature on the subject of numerical Integra- 
tions and It does not seem necessary to develop here new methods or to 
describe the old ones in detail. But it appears to be useful to 
recommend methods which have sufficient accuracy without being too 
laborious, to explain their background, and to present the necessary 
working formulas in a notation adapted to the present purposes. This 
will be done here, and for further details the reader Is referred to 
the literature. 1 

The functions which have to be integrated with respect to time 
are x^, x-^, x>>, and, if mg is not neglected, Xg. 

Let y be any one of them and assume that, for a certain time 
t = t n , y n and y n are known. A first approximation may then be 
found for the value y n +i of y at t = t n + At by assuming that 
y = y n is constant throughout the time interval: 

y n+1 = y n + ^n At (29) 


-t-See, e.g., Scarborough, James B.: Numerical Mathematical Analysis. 

The Johns Hopkins Press (Baltimore), 1930, p. 227; second ed., 1950, 
p. 2kk. 
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When this value y n+ p (in this case values of xp, xp, and x 2 ) is 
introduced in the 'differential equations, y n+ p is found, and now the 
integration may he improved hy using the average of y n and y n+ p: 

y n +i = + |(y n + ^+1) At (30) 

In figure 11 the first formula uses the shaded rectangle as increment 
and the second formula uses the shaded trapezoid, hut this result is still 
not final, because with the improved y n+ p the differential equation will 
yield another and better value for y n+ p which now should he used for the 
average. The procedure must he repeated until y n+ p and y n+ ± no longer 

change. When the time step At is well-chosen, this should occur after 
the first or second repetition. 

An example is shown in tables 3 and 4. Table 3 in identical with 
table 2 and corresponds to the simplified analysis with juq = 0. In 
table 4 there are several consecutive lines for each time t, each of 
them resulting from one complete cycle of Iteration. One sees that Xp 

and Xp are practically settled after the first cycle and that x 2 

required most of the effort. The computation may be speeded up if 
proper advantage is taken of this situation. One begins the Second line 
of table 4 with the last three columns. With x 2 = 1.287 inches and 
equations (28b) and (28a) the columns Fp, Wp - Fp, and x'p of the 
second line of table 3 may be filled. It is then possible to enter 
Xp = -155.2 inches per second squared in the second line of table 4 and to 
perform the two integrations leading to xp and xp at once with the 
trapezoid formula, equation (30). The values so obtained for x^ and 
Xp will be close to the final ones, and one may now run as many cycles 
as necessary in the Xg integration and finally check Xp and Xp 

again. Great care should be taken that the^next step is not started 
before a perfect result has been obtained, because otherwise avoidable 
errors would accumulate from step to step in the integration. 
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CD (a) 


.0025 
.0025 
.0050 
.0050 
.0075 1.879 
.0075 1.879 
.0100 2.174 

.0100 2.174 


(3) 

(4) 

(5) 

(8) 

(7) 

(8) 

(9) • 

do) 

*1 

RSI 

z 

mm 

Z 

*2 

H 

1 

Hi 

H 

\ n 

(in./sec) 


(In.) 

■ftMi 

(in. /sec) 

(in./Bec) 

(ib) 

(in./sec^) 

119.5 

0.988 

0 

12,350 

0 

119.5 

-12,350 

-119.2 

119.2 

1.287 

0 

16,090 

15.0 

104.2 

-16,090 

-155.2 

119.2 

1.288 

.018 

15,850 

14.2 

105.0 

-15,850 

-153.0 

118.8 

1.530 

.053 

19,120 

18.6 

100.2 

-19,120 

-184.7 

118.7 

1.524 

.059 

19,060 

18.4 

100.3 

-19,060 

-184.0 

118.3 

1.775 

.104 

22,200 

20.5 

97-8 

-22,200 

-214.2 

118.2 

1.772 

.107 

22,170 

20.3 

97-9 

-22,170 

-214.0 

117.7 

2.017 

.157 

25,200 

21.3 

96.4 

-25,200 

-243.4 

117.7 

2.015 

.159 

25,200 

21.3 

96.4 

-25,200 

-243.4 

• 

• 

* 

* 


• 

• 

• 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

t 

(sec) 

*1 

Aii 

*1 


mm 

nnm 

AX2 

*2 

(Iti./oec 2 ) 

(in./sec) 

(in. /sec) 

(in.) 

mam 


(in.) 

(in.) 

0 

-II9.2 

-0.298 

119.5 

0.298 

0.988 

119.5 

0.299 

0.988 

.0025 

-155-2 

-343 

119.20 

.298 

1.206 

104.2 

.280 

1.287 

.0025 

-153-0 


119.16 


1.286 

105.0 


1.268 



-.340 


.298 



.280 


.0025 


— 382 

119.16 

.297 

1.286 


.352 

1.268 

.0050 

-184.7 

-.422 

118.78 

.297 

1.583 

100.2 

.256 

1-530 

.0050 

-184.0 

-.421 

118.74 

.297 

1.583 

100.3 

.256 

1.524 

.0050 


-.460 

118.74 

.296 

1.583 


.251 

1.524 

.0075 

-214.2 

-.498 

118.28 

.296 

1.879 

97.8 

.248 

1-775 

.0075 

-214.0 

-.498 

118.24 

.296 

1.879 

97.9 

.248 

1-772 

.0075 


-535 

118.24 

.295 

1.879 


.245 

1-772 

.0100 

-243.4 

-.572 

117.70 

• 295 

2.174 

96.4 

.243 

2.017 

.0100 


. 

117.67 

. 

2.174 

• 


2.015 
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In practical computation work it is more convenient not to write 
all the lines shown in tables 3 and 4 but to erase each figure as soon 
as it can be replaced by a better one. All that is then left of table 4 
is shown in table 4(a) . As one may see, this table has the great advan- 
tage that the figures needed for averaging always stand close together. 


TABLE 4(a) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

§1 

*1 

(in. /sec 2 ) 

Ax-j_ 

(in. /sec) 

X 1 

(in. /sec) 

Ax^ 

(in.) 

1 — 1 • 

£ 

*2 

(in. /sec) 

Ax 2 
(in. ) 

X 2 

(in. ) 

D 

.0025 

.0050 

.0075 

.0100 

• 

a 

-119-2 

-153.0 

-184.0 

-214.0 

-243.4 

• 

• 

a 

-0.340 

-.421 

-.498 

-.572 

• 

• 

• 

119.5 

119.16 

118.74 

118.24 

117.67 

• 

• 

• 

0.298 

.297 

.296 

.295 

• 

• 

• 

0.988 

1.286 

1.583 

1.879 

2.174 

• 

• 

• 

119.5 

105-0 

100.3 

97-9 

96.4 

• 

• 

• 

0.280 

.256 

.248 

.243 

• 

• 

a 

O.988 

1.268 

1.524 

1.772 

2.015 

• 

a 

a 


The procedure may be accelerated considerably if at the start of a 
new step a good guess is made for the new increment instead of first 
computing a poor approximation with equation (29). If this is done and 
if the step At is chosen small enough, the method works rapidly and 
nevertheless develops good accuracy. 

Most of the criticism which this method has received in the litera- 
ture applies only to its use in problems which require a much higher 
accuracy than does the landing-impact problem. In this case slide-rule 
accuracy will always be sufficient, and this can be obtained by the 
trapezoid integration without resorting to painfully small steps. 

However, the method has the disadvantage that one never knows exactly 
how large the error is. This drawback will be avoided if the straight 
line in figure 11 is replaced by an interpolation parabola. This may be 
done as soon as four or five successive values of y have been determined. 
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Through five consecutive points (fig. 12) a parabola of the fourth degree 
may be fitted, and the coefficients of the corresponding polynomial in 
t may be written in terms of the ordinates y n _^, ^- 0 .- 3 ’ • • ■ ^n 0T} 

better, in terms of y n and of a set of differences of increasing order 


Vn = yn - y n -i 

^n = ¥n " ^n-l 


and so forth which may be computed in the following scheme: 



The polynomial may then be integrated over any one of the intervals At 
and in this way improved va l ues for the increments Ay may be obtained. 
They are computed from the following formulas: 

Ar n = - y n -i 

= At(y n - g A-j_y n - ~ 2¥ A 3 y n “ 720 ^ y n * * ’) (31a) 

Ay n _i = y n -i - y n -2 

= At ( y n - I A l y n 

Ay n -2 = y n - 2 " y n -3 

= At (y n - | A l y n 


+ 4 ^n + £7 A 3 y n + ^ Ai ^n • * •) ( 3Ib) 

23 A . 3 . _ 1 £. A . \ 

+ 12 A 2 y n _ 8 A 3 y n ' 720 * * * J 


( 31 c) 
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From these Ay's improved values of the y's are obtained. When 
they are introduced in the differential equations, better values of the 
derivatives will be found, and this procedure must be repeated until 
the results become stationary. This should occur after two or three 
cycles. If it takes longer this indicates that the time step At was 
chosen too long, and one should at once make a new start with shorter 
intervals. On the other hand, if the final .values are hit at the first 
stroke, this generally indicates that the time step was chosen too short, 
and one should start again with a greater At or continue until eight 
lines are completed and then double the step by dropping every other one. 

When this polynomial method is applied to the landing-gear problem, 
columns (1) to (6) of table 4 must be replaced by table 5, and the columns 


TABLE 5 


(1) 



( 4 ) 

( 5 ) 

(6) 

( 7 ) 

(8) 

( 9 ) 

(10) 

(11) 

(12) 

( 13 ) 

( 14 ) 

0 



Agxi 

A 

X. 

H 

A4X1 

AX]_ 

X 1 

A l*l 

aA 

A 3 ii 

Vi 

AX]_ 

x l 


- 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 


referring to x 2 by a similar table or an abridged version, depending 
on whether equations (26) or (28) are used. The results of the trapezoid 
integration are introduced into column (2), differences in columns ( 3 ) 
to (6) are computed, and then column ( 7 ) is filled with the help of 
equations (31), identifying y with x-^.. From the increments in 

col umn (7) values of X]_ in column (8) may be found which are already 
better than those of the trapezoid integration. They may at once be 
used for computing the differences in columns (9) to (12), the increments 
Axj_,- and the values x-j_, again using equations (31) • 

When this is done with the figures of table 4, it is found that 
neither the x-^s nor the X]_^s are capable of improvement, but X 2 
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is changed appreciably. The final state of the integration for x 2 is 

shown in table 6. Practically all the correction is due to the first 
time interval. 


TABLE 6 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

t 

(sec) 

*2 

(in. /sec) 

Aix 2 

A 2*2 

A^Xg 


Ax 2 

(in.) 

*2 

(in.) 

0 

.0025 

.0050 

.0075 

.0100 

• 

■ 

• 

119-5 

105.2 

100.5 

97-9 

96.5 

• 

• 

• 

- 14.3 

- 4.7 

-2.6 

-1.4 

• 

• 

• 

9-6 

2.1 

1.2 

• 

• 

• 

-7-5 

-.9 

• 

• 

• 

6.6 

• 

• 

• 

0.277 

.256 

.248 

.242 

• 

• 

0.988 

1.265 

1.521 

1.769 

2.011 

• 

• 

• 


When this polynomial method is applied - and the example demonstrates 
that it may be worth while to do so - then it will be reasonable to use 
it not only for checking and correcting but also for integrating ahead. 

To do this, one must extrapolate the polynomial in figure 12 beyond t n 
through the next interval and then integrate y from t n to t n+ p. 

The result may be expressed by the set of differences used before: 


^"n+l ^n+1 ” ^n 

/. i 5 • 251 • ' \ 

= At (y n + -A 1 y n + j^A 2 y n + § A 3 y n + — . . . ) 


( 32 ) 
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With the help of this formula one might find in table 5 the values 
of Ax-]_, x^, Ax-j_, and x^_ in the next line. But actually it is neces- 

sary to use equation ( 32 ) only once, preferably for integrating Xg 
^or Xg, if mg = 0 ), and then approximate values for the other deriva- 
tives may be procured in time to do all other integrations at once with 
equation ( 31 a) which is more exact and less influenced by the higher 
differences. 

When thus a new line in all tables (table 1 or. 2 and the integra- 
tion tables) has been filled, equation ( 31 a) is used repeatedly to 
improve x g and Xg as long as they are capable of improvement. 

As soon as the columns for the derivatives (columns (2) and ( 8 ) 
in table 5 ) fill up, one might extend the difference scheme toward 
differences of higher order, but the farther one goes to the right, the 
smaller and the more erratic the differences - will become and they will 
not be able to influence the increments computed from equations ( 31 ) 
and (32). In general the time step At should be chosen such that the 
fourth-order difference may be neglected. 

Except for the start of the computation which is always a little 
irregular, the higher differences should be rather small before they 
become erratic; otherwise one must either increase the accuracy of the 
derivatives by carrying more digits or decrease the step At. If it is 
intended to carry more significant figures, one should keep in mind that 
a many-digit machine computation is a wasted_effort, if somewhere in the 
process a figure must be read from a graph, for example, the effective 
orifice area A 3 as a function of the stroke x. 

In order to check the accuracy of the two methods - trapezoid and 
polynomial - an example of a linear shock strut has been computed with 
the following data: 


m^ = 103-6 lb sec^/in. 


mg = 0 


b = 500 lb sec/in. 

k-L = 2800 lb/ in. 
kg = 12,500 lb/ in. 


W = 0 


V = 120 in. /sec 
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The results for this example can he compared with the exact solution, 
equations (9) and (10). The following values were obtained for the 
impact force F-j_: 



F 1 

(lb) 

(sec) 

Exact 

At = 0.01 sec 

At = 0.02 sec 


Trapezoid 

Polynomial 

Trapezoid 

0.04- 

38.1 X 10 3 

38.13 X 10 3 

38.03 x 10 3 

38.^7 X 10 3 

.08 

52.0 

52.1 

52.0 

52.3 

.12 

56.2 

56.2 

56.2 

5 6 .b 

.16 

55-7 

55.8 

55-7 

55-8 

.20 

52.5 

52.6 

52.6 

52.6 

. 2 b 

bj .6 


bj .6 



Evidently, under these conditions the trapezoid method with 
At = 0.01 second is good enough. Encouraged by this result, the 
step At has been doubled. The results of the trapezoid integration 
are shown in the table. The polynomial method proved to be extremely 
tedious and was not pursued further when after several hours of computa- 
tion the first four lines had not yet stabilized. However, it was found 
practical to start with the small interval and double the step as soon 
as possible. The polynomial method with At = 0.01. second was carried 
to t = 0.08 second, and then the results for t = 0.02, 0.0^> 0.06, and 
0.08 second were used to start the polynomial method with the double 
interval. This computation waB carried up to t = 0.20 second and 
yielded results identical with those obtained for the shorter steps. 


Numerical Example 

As an illustration of the methods just described an example has 
been worked out. The data chosen and the metering pin correspond closely 
to those of a recent American airplane. The data are these: 


A-j_ = Ag = 39-8 sq in. 
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Vo 



= 23-5 in. 


p Q = 310 Ib/sq in. 
p = 8A2 X 10“ 5 lb sec 2 /in.^ 


7 = 1.1 


kg = 12,500 Xb/in. 


W-l = 0 

V = 10 ft/sec 


The effective orifice area is shown by the heavy line in figure 13 
as a function of the stroke x. The low part at the left-hand Bide of 
the diagram represents the bulbous end of the metering pin. 

The shock strut is prestressed with the. force 


F 0 = p Q Ag = 12,350 lb 

Until the impact force has reached this value, only the tire is deformed 
and the simple formulas mentioned after equations (2 6) apply. They 
yield t 0 = 0.0088 second, X 0 = O.988 inch, and V 0 = 119-5 inches per 
second. These are the initial conditions for the numerical integration 
of equations (25) and (28). This integration was started by the 
trapezoid method, using equations (29) and (30), and the time step At 
was so choseh that at least a few intervals would pass before the first 
break in the curve A3 = A^(x) was reached. This is possible with 

At = 0.0025 second, and the first lines of this computation are shown 
in tables 3 and 4 (where t is counted from the beginning of this Integra 
tion, not from the first contact between tire and runway) . When four 
steps were completed, the polynomial method was started and the results 
of these steps were improved. The computation was carried on to 
t = 0.0250 second with x = 0.571 inch, X]_ = 3*907 inches, and 
Xg = 3.336 inches. . This is sufficiently far past the first break in 
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A^(x) that it was possible to double the step. A new integration table 

was started with the results for t = 0 . 010 , 0.015, 0 . 020 , and 0.025 second 
and it was carried on with At = 0.005 second until t = 0.060 second with 
x = 1.939 inches. The next step would have led beyond the second break 
in the curve A 3 (x) and hence to large values in the difference schemes. 

Therefore, it was necessary to return to the shorter time step 
At = 0.0025 second. This makes it necessary to interpolate values for 
the half intervals. To keep up with the accuracy of the integration, 
this must be done with the help of the same interpolation parabolas from 
which equations (31) and (32) are derived. With the notations used the 
following formula holds: 

y n+ 1 = |(y n + y n +i) - jgfe & n+i + ^aw) + ^(^+2 + ^+3) (33) 


in which the last term is often negligibly small. 

With the help of this formula a new integration table was started, 
beginning with t = 0. 0^50, 0.0475* . . . second. When it came to 
t = 0.0775 second, the next and last break in the A 3 diagram was 

reached and the higher differences rose so high that it became necessary 
to reduce the step to 0.00125 second. Eight lines beyond the discon- 
tinuity the step was increased to At = 0.0025 second and soon thereafter 
to 0.005 second. At t = O.lh second it was realized that the higher 
differences had become so small that the interval could again be doubled, 
and with At = 0.01 second the computation was carried until 
t = 0.27 second, when x became negative. 

The example which was chosen here as a test specimen for the 
numerical integration is one of the most irregular possible. Most of 
the computation effort was spent on the bulbous end of the metering pin. 

As soon as the last corner in figure 13 was passed, the work proceeded 
rather quickly to its end. When the pin is shaped more gently, or when 
there is no pin at all, it will be possible to start, say, with 
At = 0.005 second and to change after some time to At = 0.01 second, 
without the many tedious changes which -were necessary in the present case. 

The results of the computation are shown in figures lb and 15 • There 
is a double time scale in the diagrams, one beginning at the first contact 
and one at the time t 0 when the numerical integration begins. 

Figure l4 shows the stroke x and the displacement of the 

airplane. There is a first, short phase during which only the tire is 
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deformed and x = 0 . Then the shock strut begins to work but, because 
of the bulbous end of the metering pin, the strut collapses, at first 
rather slowly . Later it catches • up, and the curves x and xj_ approach 
each other, indicating that the load maximum is passed and that the tire 
expands . 

Figure 15 shows the load and its breakdown into the damping force 
and the elastic (air) force according to the two terms of equation (25) • 
Because of the bulbous end of the metering pin the damping force builds 
up rapidly, but then the orifice opens up and the increasing air pres- 
sure in the shock strut cannot compensate the decline of the damping 
force. 


Dimensioning of Metering Pin 

For reasons of weight saving it is desirable that the shock-strut 
force rise quickly to a high value and then remain at this height for a 
sufficient time to bring the mass m^ to rest. As a practical means 

for this purpose, the metering pin has been introduced into the design 
of shock struts. Wow, since there is but one metering pin, it will not 
be possible to obtain ideal results for different impact conditions, but 
it is possible to pick out one landing case of particular importance and 
to shape the metering pin so that in this case a desired load history is 
obtained. The shape of the pin which has been found for this case must, 
of course, be subjected to a critical study in two respects: It must be 

acceptable to the workshop, and it must yield at least tolerable load- 
time diagrams under other landing conditions. Th* final compromise is. 
a true engineering decision which cannot be replaced by an analytical 
device. 

There is no need to specify exactly how the impact force should rise 
from zero up to a certain level. In this first part of the load history 
the tire has an important influence, and it will be enough to choose the 
orifice opening A3 so that not too much stroke is lost while the force 
builds up. 

But when at a certain time t = t' the force F^ has reached a 
certain value, say Fj_ =-F', then it may be desirable to keep it constant 
on this level. If it is agreed to neglect the unsprung mass m£ equa- 
tions (28) are simply 


m^xj + F' = Wj_ 


(3^a) 

(34h) 


F' = kgXg 
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and it follows at once from equation (3^£>) that Xg = Constant, say, 

Xg = Xg' 

Xg = 0 

Equation (3^a) presents a simple integration problem and yields 

W, - F' 

Xp = Xp ’ + ( t - t * ) 

m l 


= Xp' + Xp'(t - t f ) + 


W, - F' 


2nu 


-(t - t *) 1 


where ip* and xp' are the values which the variables have assumed 
at t = t'. From equation (28c) 


x = 


X = Xp - Xg' 


and these values may now be introduced in equation ( 25 ) of the oleo strut, 
which then yields 



A lP c 


/ z A ' 
I O 0 

(z 0 A o - xAl 


7 


(35) 


This idea has been applied in two ways to the numerical example of 
the preceding section. 

When looking at figure 15, one might think it useful to keep Fp 
for some time on its peak level, thus decelerating faster the vertical 
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motion of the airplane without imposing a higher dynamic load on it. 

This would result in a saving in stroke and "hence in weight of the 
shock strut. In this way the curves marked "I" in figures 1 6 and 17 
have been obtained. The corresponding orifice area is shown in fig- 
ure 13 by “the line. I. To keep the impact force at its peak level, the 
orifice area must be decreased relative to the original design, and it 
comes down to zero, when the motion of the airplane is stopped. The 
steep descent at the end of this curve is, of course, not acceptable 
for the design, since it. means a complete - plugging of the orifice and 
would lead to a high load peak in a case ofyharder landing) but the 
upper part of the curve may lead to an improvement of the design. 

One might think of another modification., of the load-time curve, 
cutting away the peak and. keeping F]_ as long as feasible on a medium 
level, say at F* = 46,600 pounds. When this' is done, the curves 
marked "II" in figures 13, l6, and 17 result. They show that in this 
case a slightly longer stroke is needed than in the original design, 
but there is a considerable saving in dynamic load. 

Since the rise of Fj_ is interrupted in this case, the orifice 
must be opened wider, and figure 13 shows that most of the bulbous end 
of the pin must be removed. The transition must, of course, be smoother 
than that shown in the diagram, and this would lead to a rounding of the 
corner in the load-time diagram (fig. 17)- Except for this necessary 
modification and for the steep end of the A3 curve, the solution seems 
acceptable, provided that the pin shaped in this way proves to be satis- 
factory in other landing cases. - 

But there is still one essential point that needs discussion. 

Figure 18 shows the velocities Xj_ and x for all three cases. For 

the original pin heavier lines have been used and the two modifications 
are marked "I" and "II." The first modification does not show anything 
in particular, but for the second modified pin x jumps suddenly from 
one value to another and so does = x^ - x. Now, a sudden change of 

the velocity X£ will, of course, meet with the inertia of the unsprung 

mass, and the metering pin II cannot be accepted without discussing this 
influence. " — - ■ 

Starting from equations (26) and putting F^ = F ' : 


mixi = W x - F' 


mgXg + kgXg = F'„ 
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The first of these equations is identical with equation (3^a)j and the 
second yields an undamped vibration: 


Xg = + A cos os(t - t 1 ) + B sin o>(t - t') 


03 = 



Now, at t = t f , when this vibration begins. 


Xg = Xg’ = F'/kg 



x 2 


and hence 


Xg = Xg’ 



sin cu(t - t’) 


When it is assumed in the example that mg/m p = 0.025, the circular 
frequency of these vibrations is o> = 69-5 second -1 , that is, about 
10 cycles per second. 

The stroke x will show the same undulation as Xg and so will 

the metering pin. Of course, nobody would think of building a metering 
pin of that shape, in particular since the length and location of these 
undulations would depend on the arbitrary choice of the conditions under 
which Fp is kept constant. As soon as a streamlined metering pin is 
chosen corresponding to the simplified analysis, the force F^_ will 
fluctuate slightly and thus provide the necessary damping for the transient 
vibrations of xg and x. 

There is still a better way of handling this last question. Since 
it is not feasible anyway to make a metering pin with a sudden change of 
cross section, it is better to assume a force diagram on which the corner 
is well-rounded, say by a parabola 

p 

F^ = c^ + Cgt + c^t 

which is so chosen that there is no large discontinuity in dF^_/dt. 
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When this force-time relation is introduced into the equations of 
the landing gear, they may easily he integrated, and the resulting 
expressions for x and & along with F may he introduced into equa- 
tion (35) to find A3 and hence the cross section of the metering pin. 

The short broken line at the outset of .the horizontal.. line marked 
"II" in figure 17 represents such a parabolic rounding of a corner in 
the force diagram. The corresponding values of x and A3 have been 

indicated by broken lines in figures 18 and 13 • One may recognize that 
no great change of the metering pin is needed to make the wheel motion 
much smoother, and the corner in the force diagram might still be rounded 
much more without a substantial loss of deceleration for the airplane. 


ADDITIONAL NONLINEARITIES 
Tire 


The elastic resistance of the tire depends only in small part on 
the elasticity of the rubber and is essentially due to the compression 
of the enclosed air. During the landing impact this compression is 
nearly adiabatic and therefore the relation between the tire pressure 
and the deflection x 2 is nonlinear. On the other hand, the relation 

between the pressure and the force F 2 is nonlinear also because the 

tire flattens. On the whole, these and some other influences seem to 
compensate to some extent, and load-deflection curves from tests may be 
fairly well approximated by a straight line.- -This is illustrated by 
figure 19 which shows such a test result. 

For design purposes it does not seem worth while to replace, under 
these circumstances, the linear relation (equation (la)) by anything more 
complicated. However, for the. evaluation ofutests it may be advisable 
to use the best available information on the behavior of the tire. 

The nonlinearity of the tire becomes severe when it comes to 
bottoming. . Then the force F 2 may rise to high values without an 

appreciable further increase. of x 2 - In general, bottoming should, of 
course, be avoided, but when it comes into .consideration, then equa- 
tion (la) can no longer be applied, and it must be replaced by the general 
relation ~ 


F 2 = F 2 (x 2 ) 
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which represents an empirical function determined from tests. The 
equations of motion are then these: 


+ ^i(x, x) = 


(36a) 


nigXg “ ®i(x, x) + Fg^Xg^ = 0 


(36b) 


instead of equations (26a) and (26b). 

Because of the prestressing of the oleo strut these equations are 
not valid until the shock-strut force has reached the prestress value F Q . 
For this initial phase of the impact the procedure described in the 
paragraph following equation ( 26 c) must be applied. Since it covers 
but a small part of the whole impact, one may use there the linear 
law F2 = k2X2, the spring constant kg being taken from the initial 
tangent of the load-deflection curve of the tire: 


kg = 


^2(^2) 

dXg 


x 2 =0 


For equations (36a) and (36b) then the initial conditions (equa- 
tions ( 27 )) are the same as those for equations (26). The equations are 
solved by numerical integration and table 1 may be adopted, changing 
only the heading of column ( 9 ) where Fg is written instead of kgXg 

and then using a graph of the function Fg^Xg) to fill this column. 

In most cases it will be possible to neglect the mass mg. Then 
equations (36) are rewritten in the form 

m ix*i = W x - Fg^Xgj ( 37 a) 

F x (x, x) = Fg(xg) ( 37 b) 


which corresponds to 
table 2 , writing Fg 
column with the help 


equations (28). For the numerical integration use 
= F^_ at the top of column (6) and filling this 
of the graph for Fgfxg). 
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In both cases, equations ( 36 ) and (37)> the integration step must' 
be decreased appropriately when approaching; the region where the tire 
bottoms. 


Kinematic Nonlinearities 

In the equations of motion one needs the second derivatives of the 
displacements x^_ and X 2 , the accelerations of the masses m^ and 

respectively. For the shock strut, the stroke x is needed. Thus far 
it has always been assumed that x is equal to the difference x^ - Xg. 

However, this relation holds only in the simple case, when the upper part 
of the shock strut (usually the barrel) is rigidly connected with the 
airframe and the wheel is attached directly to the lower part (piston). 

A correction is already needed when the shock strut is inclined from 
the vertical (fig. 20). In this case 


x i - Xg 

x = — - 

cos a 

The changes which this relation requires in the intergration schemes 
are obvious and there is no need to discuss them in detail. 

However, there are cases in which the relation between x and 
(xj_ - Xg) is nonlinear. Figure 21 illustrates what is meant. Most of 

these devices have disappeared from current -practice, but in a time of 
rapid development it is advisable to discuss briefly how similar cases 
may be handled. For all these landing gears a nonlinear relation 


X = f( X] _ - x 2 ) 


( 38 ) 


can be established by trigonometric methods. By differentiating it with 
respect to time, the relation 


. df (*i - 

d(xi - Xgj\ 1 

= f'( Xl - Xg^x-L 




(39) 


is derived. The above two equations take the place of equation ( 26 c) 
and the corresponding relation for the velocities. 
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Table 1 must now be replaced by table J } in which, columns ( 5 a), 
(6a), and (6b) have been added. Columns (5a) and (6b) are filled from 
the preceding ones and column (6a) is filled from a formula or a graph 


table 7 


(1) 

(2) 

( 3 ) 

w 

( 5 ) 

( 5 a) 

(6) 

(6a) 

(6b) 

( 7 ) 

(8) 

( 9 ) 

(10) 

(11) 

(12) 

t 

X 1 

a 

*2 

*2 

X 1 “ *2 

X 

f* 

*1 ‘ *2 

X 

F 1 

Va 

H 

w i - F i 

*1 

















for f'(x-L - x 2 ). Then equations (38) and ( 39 ) are used to fill 

columns (6) and (l) . Everything else is done as explained for 
table 1. 

When m2 is neglected, table 8 is used instead of table 2 . Again 
additional columns for ^x^ - Xg), (x-^ - Xg), and f* are provided. 
Column ( 5 ) is filled with the help of equation ( 38 ); column ( 7 a), from} 


TABLE 8 


(1) 

(2) 

( 3 ) 

( 4 ) 

(*«) 

( 5 ) 

(6) 

( 7 ) 

( 7 a) 

( 7 b) 

(8) 

( 9 ) 

do) 

0 

X 1 

X 1 

*2 

X 1 ' *2 

X 

F 1 

X 

fl 

X 1 - X 2 

*2 

9Q 

' x i 















column (ka) with a graph or formula for f'(x]_ - x 2 ); and column (7b), 
with equation (39)- 

These are very simple changes, the numerical integration being a 
very flexible instrument that can be adapted to almost every special 
requirement . 




























6o 


NACA TN 2743 


THE AIRPLANE AS A WHOLE 
Introduction ~ 


In most investigations of landing gears the airplane is represented 
by a single mass m-j. aiding on a system of .springs and dampers, with 

perhaps a small additional mass mg representing the wheel. All the 

preceding sections of this report are exclusively concerned with this 
model. 

However, the real airplane is a three-dimensional structure, and when 
one or more of its wheels hit the ground, it may receive not only a 
vertical acceleration but also angular accelerations about different 
axes. These angular accelerations and the rotatory motion resulting 
from them will, of course, influence the landing impact. 

A detailed study of this phenomenon leads into rather lengthy 
computations . Their quantitative results will depend on many details 
and may vary widely between different types of airplanes. This section 
will therefore be restricted to some general^considerations concerning 
the best method of analysis. 

There are two principal problems, the symmetric case in which both 
wheels of the main landing gear strike the ground simultaneously and in 
identical conditions, either earlier or later than the auxiliary gear, 
and the asymmetric case in which the two wheels of the main gear touch 
the ground one after the other. 


Symmetric Impact 

Figure 22 shows the side view of an airplane as far as it is of 
interest for the present purposes. The point - C is the center of gravity 
where the mass m is located. To the right is the main gear; to the 
left, the auxiliary gear which may be either a nose gear or a bail gear. 

In figure 22 the airplane is shown In the position which it has at 
the time t - 0, when the main gear makes its first contact with the 
runway. From this time an impact force F^ of increasing magnitude will 

act in each main gear and it will cause both’ a deceleration of the 
vertical movement of the center of gravity and a pitching motion about 
this point. The equations of motion are 
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mz = -2F ± + W 
Iy(b = -2F 1 a 



> 


j 


(to) 


where I is the moment of inertia of the airplane with respect to its 
transverse axis. 

The resultant acceleration at the upper end of the main landing 
gear will he 


x-, = z + am 


f 1 ¥ 

m - Iy m 


EiA + sd) + « 

m \ Iy I m 


In the section "Linear Spring -Damper Systems" there was written 


F 1 + W 1 
Xt = - + 

i m^ 


and the two expressions are equivalent if one chooses 


m, = 


m 


2 



V 


m 

2 2 .2 

a +1 y 
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where 



is the radius of gyration. 


These formulas show that 


it is perfectly justifiable to study an isolated landing gear, provided 
one does not simply use as mass m-j_ one-half of the airplane maBs. 

However, this procedure is subject to two essential limitations: It 

can be applied only until the auxiliary gear comes into action and must 
at least be modified when the rotation of the airplane leads to a sub- 
stantial change of the angle of incidence of the wing and hence to a 
change of the load W. 


Consider the second point first. The angular position of the air- 
plane is determined by the angle a between the ground and a reference 
line in the plane of symmetry of the airplane. This reference line is 
so chosen that a = 0 when all three wheels - of the airplane just touch 
the ground without pressure. The angle a which is so defined is not 
identical with the angle of incidence of the wings, but the two differ 
only by a constant which depends on the design of the airplane. 

Since only small values of a need be considered, it may be 
assumed that the lift and hence W is a linear function of a, say: 


W = W' + W"a 


but since one must use numerical integration methods anyway an arbitrary 
function 

W = W(a) 

may be assumed when this appears to be necessary. The part of this 
weight which must be attributed to one main landing gear is then 


Wi(a) 



( 41 ) 



NACA TN 27^3 


63 


Equations of motion (40) are now written in the following form: 


*1 - - 5^*1 - w i) 


(42a) 



(42b) 


Additionally, there is a relation which connects F-, with x 

and x, for example, equation (25) of the oleo strut, the relation 
x = Xj - xqj and the elastic equation of the tire Eg = Fp = 1^X2 • 

Table 2 which is used for the one-gear problem must now be extended 
so that it may take care of equation (42b). It looks then as shown 
in table 9: 


TABLE 9 


(1) 

(2) 

(3) 

w 

(5) 

(6) 

(7) 

(8) 

(9) 

( 10 ) 

(ID 

(12) 

(13) 

(14) 

t 

X 1 

X 1 

*2 

X 

a 

• 

a 

*1 

• 

X 

*2 

W 1 

W 1 -Fi 

• • 

X 1 

• • 

a 

0 













• 


Columns (2) to (5) and (8) to (10) are treated exactly as are the 
corresponding columns of table 2; the first line in co lumn s (6) and (7) 
is filled in from two additional boundary conditions (a given, a = 0) . 
Also the starting value of Wi will be known and must, of course, check 
with column (6) and equation (4l) . Column (12) is self-explanatory, and 
columns (13) and (l4) follow from equations (42). Besides the tables 
for the integration of Xp, Xp, and Xg, an additional table is now 
needed to integrate a and a. With the results of these integrations 
the second line may be started. 

Of course, this analysis does not consider the possibility that the 
pilot uses the controls to counteract the pitching movement of the air- 
plane. If he does bo, a human element comes into play which is not 
easily incorporated in mathematical formulas. This uncertainty may 
upset the usefulness of the procedure and will justify the application 
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of the simpler table given previously. This simplification, commonly 
used in landing-gear analysis, is still more justified by the results 
represented in figure 6 , which show that the exact magnitude of the 
effective weight Wp is of secondary importance for the interesting 
portion of the impact. 

Whether a simple one-gear analysis is made or the variability of 
W]_ is taken into account, this computation ends at the moment when the 
auxiliary gear comes into action. The time t = t' when this occurs 
is found in the following way: During the first phase of motion the 

acceleration at the upper end of the auxiliary gear is: 



When it is assumed that the airplane approaches the ground with the 
vertical velocity V, but without an angular .velocity, then the velocity 
x^ for t > 0 is 



Under the integral signs and W must be introduced as functions 

of - t according to the analysis of the main gear. 

The displacement x^ of the auxiliary gear is best counted from 

the position in -which the wheel just touches the ground. When the 
airplane lands at an angle a (fig. 22 ), then x^ = -(a + b)a at 

t = 0. For t > 0. 

"t . ■ 

Xo = ~(a + b)<x + x dt (44) 

Jo J 

and the time when this equation yields X 3 =0 is the time t = t'. 
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p 

When W = 0 and ah < 3y , the velocity x^ will decrease and 

the impact of the auxiliary gear will he softer than it would he if 
this gear had hit the ground before the main gear. If W ^ 0, there 
is an additional positive term in equation (43) and, since the airplane 
is still falling under the influence of the force W, the velocity X3 

may increase. When ah > iy^, then x^ will certainly increase, pos- 
sibly even very much so, and the auxiliary gear may strike the runway 
rather forcibly. 

When all wheels are In contact with the ground, the equations of 
motion are rather involved. When all landing gears have a spring- 
damper unit as a shock strut, no damping in the tire, and no unsprung 
mass, the problem is of sixth order. It is of little value to establish 
the formulas for the linear case, but it is useful to develop a numerical 
procedure which may be applied in linear as well as in nonlinear cases. 

The equations of motion contain now the forces in main and 
auxiliary gears (fig. 23): 


mz = -2Fj|_ - + W 

► 

IyOi = -2F-j_a + F^b 

Then there are two kin ematic relations 


(45) 


X 1 
x 3 

Differentiating there and then introducing z and <h 
yield 



= z + aoj 


= z - bco 


(46) 


from equations (45) 


( 47 ) 



66 


NACA TN 2743 


These equations may he used in the following way: For each landing 

gear a single-gear analysis is started according to the instructions 
given in the section "Intentional Nonlinearities . " For the main gear 
it begins at t = 0 and runs exactly as explained there until t = t'. 

For the auxiliary gear it begins at t = t 1 _with x^ = 0 and the 

value of x^ which follows from equation (43). When it is not desired 

to neglect the unsprung mass, table 1 is used, otherwise table 2. In 

either case the line for t = t' may be filled up to the last two 

columns, but the last two columns are replaced by some col umn s which are 
adapted to equations (47). They yield, x^ and x^ in terms of the 

forces F of both tables, and these values are now integrated just as 
was done with x-^ in table 4. 

A step-by-step integration of this kind requires twice as much 
time as a single-gear analysis and will yieid everything needed for 
both gears. 


One -Wheel Landing 

It is possible that a landing airplane may approach the runway 
with one wing low and that the wheels of the main gear do not hit the 
runway at the same time (fig. 24). There are then again two phases, a 
first one while only one wheel is in contact with the ground and a 
second one when both wheels are. 

In the first phase there is only one force F^, having the 

distances a and c from the lateral and longitudinal axes, respec- 
tively. It produces the following accelerations: 

Vertical at center of gravity: 


mz = -F x + W 


Angular with respect to lateral axis: 

Vv = - F i a 


Angular with respect to longitudinal axis: 
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The resulting acceleration at the upper end of the active landing gear 
is 

x-|_ = z + auy + cdijj. 

and when use is made of the preceding equations there is obtained 



Again it is useful to introduce the radii of gyration by: 

U 2 - V m 

V 2 ■ V" 

and to write 



When this is compared with the relation 



used previously, it is seen that one must put 

m 
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With these notations the' impact problem is again reduced to a one-gear 
problem until the other wheel meets the ground. The time t = t" at 
which this will occur may be found in the same way as in the case of a 
two -point landing. 

The acceleration at the top of the second main gear is 


x*2 = z + any - coy 



Integrating once yields the velocity: 




rt 

w 

F-, dt + - t 
x m 

Jo 


and integrating again, 

^ . 

x -3 = -2c|3 + Xo dt 
J JO 2 


At the time t = t" when x^ = 0, the one-gear problem ends and 
from then on both main gears must be dealt with simultaneously. 

This is done as in the preceding section, but the formulas differ 
in details because there is still one degree of freedom left, the 
rotation about the transverse axis of the airplane. During this phase 
of the landing impact the equations of motion (fig. 25) are as follows: 


m z = -Fj - F^ + W 


VV - -(*1 + F 3) a 

I A - (-*1 + f 3 ) c 
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Besides there are the kinematic relations 


Xj = z + adiy. + cd^ 
x '2 = z + adiy - cdijj. 

and by combining both sets of equations the following equations are 
obtained which correspond to equations (47): 



> m 



These equations may be handled exactly in the same way as equations (47) , 
with, however, the restriction that the auxiliary gear must still be 
off the ground. As soon as it makes contact, the relations become more 
involved, but it seems at present not necessary to elaborate the details 
of the third phase of the impact which then will follow. 


REVIEW OF GERMAN LITERATURE ON LANDING-GEAR IMPACT 


Before the last war in Germany almost no theoretical work was done 
on landing-gear problems, and it seems also that in other countries 
interest was low. 

During the war in Germany new and unexpected demands could fre- 
quently best be met by adapting an existing airplane type, with its 
well-established mass -production facilities. Such modifications usually 
resulted in an increase of weight without supplying additional space 
into which a larger wheel could be retracted. Frequent tire troubles 
were the unavoidable consequence, resulting in a strong impetus to 
landing-gear research. All but one of the papers reviewed here belong 
to this period of wantime research. 

When studying this German wartime literature, one must keep in 
mind during what period and under what circumstances the work was done. 
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All of these papers appeared in 1943 and 1944 and were thus the outcome 
of a rather short period of research. They represent an intensive 
attempt to tackle a" long-neglected problem. But before the work had _ 
yielded results of “final validity, it was cut off early in 1945 by 
the national catastrophe. More than 6 years have elapsed and the 
landing-load problem has undergone changes. Some of the statements 
made in those papers have lost interest, others are no longer applicable 
without modification, and most of the analytical methods are either 
oversimplified or too complicated. 

Nevertheless, it is still worth while to survey this literature 
briefly because it “contains 'many of the ideas and methods which are 
still the basis of . landing-gear analysis. Indeed, in writing this 
report the author has drawn much useful information from the German 
publications which^are reviewed on the following pages. 

The goal of the early landing-gear research was influenced by the 
attitude of official regulations. They required that a drop test be 
made in which the upper end of the shock strut was connected with a 
mass in the notation of this paper) and the two dropped on an 

anvil. At the instant when the anvil was struck, the weight Wp was 
compensated by admitting compressed air to two cylinders. The load- 
stroke curve obtained by this test was then considered as "the” load- 
stroke curve of the shock strut and was employed in all landing cases 
which had to be considered in the design of the airplane. _ Consequently, 
the effort x>f the early research was directed toward the investigation 
of load-stroke diagrams of shock struts. II 7 

The first paper that must be mentioned here, and the only one that 
appeared before the war, was written by Michael (reference l) . It gives 
a detailed analysis of ..the linear spring -damper system but pays only 
slight attention to the tire. A special feature of this paper is the 
use of spring diagrams in which the f.orce is plotted either against the 
stroke with the rate of stroke as a parameter .or, inversely, against 
the rate of stroke with the stroke as parameter. These diagrams are 
shown also for shock struts with dry friction or with velocity-square 
dampers, and they are used for a graphical solution of the differential 
equation. Such diagrams are no longer possible when a second spring 
(the tire) is present, and therefore they have not been employed again 
in later papers. 

The first papers of the war period were .still focused on the load- 
stroke diagram. Schlaefke (reference 2) criticized the drop-test method 
and suggested replacing the buffered drop test by an unbuffered test, 
that is, omitting the air cylinders and with them a possible source of 
inaccuracy. His paper uses the theory of the linear spring-damper system 
to establish some relations between the results of both tests. 
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In a later paper (reference 3) the same author realizes that the 
.damping in the oleo strut is far from proportional to the rate of 
stroke. He compares load-stroke curves for lin ear and for velocity- 
square damping and arrives at the strange conclusion that the former 
look more realistic. The method used for the analysis of the nonlinear 
problem is of interest. A balance of kinetic and potential energy is 
established and from it, a differential equation between x as dependent 
and x as independent variable. When it is solved, the damping force 
(proportional to x2) is also known in terms of x. However, this 
approach is not possible in the presence of a tire. 

In the next group of papers the tire makes its appearance. A 
paper by Kochanowsky (reference 1+) gives a very detailed analysis of 
the oleo-tire combination as shown in figure 1. Kochanowsky finds that 
the unsprung mass is of no great importance for the landing impact nnd 
that the problem may readily be simplified by assuming mg = 0. The 
study of this paper (and of many others) is rendered difficult by the 
author's habit of using for all and everything dimensionless quantities 
so that the reader has to learn first a system of not very suggestive 
notations before he can follow the analysis or read the diagrams. 

Another paper by Schlaefke (reference 5) covers approximately the 
same ground. 

After having studied the linear oleo-tire system, the next logical 
step would have been to consider a nonlinear shock strut, but, inciden- 
tally, the few papers which did this were older than Kochanowsky 's com- 
prehensive paper on the linear system. One of them is by the same author 
(reference 6), and it was not thought to be a study of a nonlinear case. 
It is concerned with a special type of spring which haB long been used 
in railroad-car bumpers and was introduced in landing gears. It con- 
sists of a pile of rings with conical sides (fig. 26). When it is 
subjected to an axial compression, the hoop stresses in the rings are 
alternatively tensile and compressive. During the elastic deformation, 
the rings slip on one another and the pile becomes shorter. Because 
of the slip, there is considerable dry friction, and when the load F 
is decreased, the deformation x is not immediately decreased but 
follows a law which is described by figure 27- The area of the tri- 
angular loop represents a loss of energy and the ring pile may thus be 
used as a damped spring. Kochanowsky 1 s paper considers a shock-strut 
and tire combination in which the strut has no other elastic or damping 
element except such a ring -pile spring. During the first upstroke the 
analysis is extremely simple, since not even damping appears explicitly 
in the equations^ but when the motion is followed beyond the force 
maximum, it is linear only in sections but nonlinear on the whole. The 
paper is an interesting study, but the ring-pile shock strut is not 
versatile enough to stand the competition with the modern oleo strut, 
and the problem is now obsolete. 
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The other paper which considers a nonlinear shock strut is a very 
serious and very detailed study by Marquard and Meyer zur Capellen 
(reference 7 ) . The authors consider velocity-square damping and 
polytropic compression of- the air, formulate differential equations, 
and integrate them numerically. Unfortunately, the authors overestimate 
the accuracy requirements of the analysis. -In their tables values are 
given to six and even seven significant digits, and consequently they 
employ an exact but very tedious method of .step-by-step integration. 

In addition to the detailed treatment of the nonlinear shock strut, 
the paper is remarkable' for another reason. It not only considers an 
oleo-tire system with a very realistic shock strut, but it also considers 
the motion of the whole airplane in its plane of symmetry. In a second 
paper (reference 8 ) the same authors extend_their investigation to cases 
of unsymmetric landing. But here also the attempt at exactness goes too 
far when the decrease of horizontal speed during the short impact time 
is taken into consideration. This is pointed out in a paper by Schmitz 
(reference 9). This author also considers the pitching motion of the 
airplane and includes the ensuing change of the lift, but he falls back 
to the old idea of "the" load-stroke curve and fails to realize that 
the cooperation between the elastic reaction of the air and the damping 
force caused by the orifice depends largely _on the conditions of the 
impact . 

Besides the landing impact, the taxying of the airplane has always 
met with interest. Michael's paper (reference 1) pays attention to it, 
and Kochanowsky ' s papers (references k and 6 ) both consider the taxying 
impact in full detail. In these papers the_ statement is made and proved 
that when the airplane rolls over a sinusoidal ground swell, the mass m-^ 
travels practically on a level path and that. therefore the analysis may 
be made on the- as sumption that mp = «>. 

Besides these papers there are two by Schlaefke In which taxying 
is considered. One of them (reference 10) covers the same ground as 
the corresponding part of Kochanowsky 1 s paper (reference 4). The second 
(reference ll) Is a short note concerning the impact during the take-off 
r un . It seems to. be the only paper devoted, to this subject, and not 
much information is found in it. 

Additionally, there are a number of reports on experiments. .Most 
of them were tests made by the airplane manufacturers and served 
essentially the purpose. of improving a new airplane model to the point 
where it was ready for production. Today, it is difficult, if. not 
impossible, to draw other than qualitative ^information from these reports 
since the- airplanes, shock struts, and tires used in these tests no 
longer exist and details needed for an analysis may no longer be obtained 
readily. 
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However, one of the experimental papers must be mentioned in this 
review, a short report of Hoke (reference 12) on the experiments he 
made in the Deutsche Versuchsanstalt fur Luftfahrt. He measured, as 
functions of time, the vertical velocity of the airplane immediately 
before and during the landing impact and the vertical and lateral forces 
on the wheel. The fine experimental technique of the velocity measure- 
ment is described in the paper. 


Stanford University 

Stanford, Calif., November 15, 1951 
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Figure 2.- Mechanical system in two positions. 
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Figure 3*- Influence of unsprung mass. 
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Figure 5-- Concluded. 
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Figure 13 •- Cross section of metering pin. 
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Figure iA.- Displacement Xj and stroke x against time. 



Figure 15 •- Force-time history of shock strut. 
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Figure 18.- Vertical velocity x-^ and rate of stroke x for different 

metering pins. 



Figure 19 •- Load-deflection curve of a tire. 
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Figure 24.- Schematic views of airplane during one-wheel landing. 

Situation at first contact. 



Figure 25-- Schematic views of airplane during one-wheel landing. 
Situation when second wheel hits runway. 




